NON-FUCHSIAN EXTENSION TO THE PAINLEVE TEST 



Micheline Musette and Robert Contef 

Dienst Theoretische Natuurkunde, Vrije Universiteit Brussel, 
B-1050 Brussel, Belgique 

f CEA, Service de physique de l'etat condense, Centre d'etudes de Saclay, 
F-91190 Gif-sur-Yvette Cedex, France 

Short title: Non-Fuchsian Painleve test 



PACS: 0230, 0340 K 



Keywords 

singularity analysis 
Painleve property 

Abstract - We consider meromorphic particular solutions of nonlinear ordi- 
nary differential equations and perform a perturbation a la Poincare making their 
linearized equation non-Fuchsian at the movable pole and Fuchsian at infinity. When 
the nonlinear equation possesses movable logarithms, they are detected sooner than 
with the perturbative (Fuchsian) Painleve test. 



Phys. Lett. A, received 24 February 1995, revised 1 August 1995, accepted 4 August 1995 

Saclay SPEC 94/118 



1 



1 Introduction 



Nonlinear differential equations with fixed critical points define a natural exten- 
sion |l] of linear equations. Let us recall that a singular point is said critical if several 
determinations of the solution are permuted around it, and movable (contrary fixed) 
if its location depends on the initial conditions. The Painleve property (PP) of a 
differential equation, is defined [l|] as the absence of movable critical points in the 
general solution of the differential equation. 

The Painleve test || is the set of all methods to build necessary conditions for 
a differential equation to possess the PP, without guarantee on their sufficiency. 
The most powerful of them is the "a— method" of Painleve j|, but its difficulty is 
differential since at each step one must integrate a differential equation. All other 
methods (Kowalevskaya [|, Bureau 0, Conte, Fordy and Pickering |7|, [J) have 
only an algebraic difficulty and therefore are easy to automatize. 

But all those algebraic methods share a restriction, which this Letter proposes to 
remove, thus adding to the Painleve test a new algebraic method which very often 
allows to conclude to a failure more rapidly than with the old algebraic methods. 

In section ^|, we recall the lemma of Painleve which establishes the perturba- 
tive framework, then we define our extension to the Painleve test. Section |^ is a 
pedagogical example quite useful to explain the method. Two examples are then 
considered, one from mathematics section f|, the second one from physics section 
Both exhibit movable logarithms and provide a shorter proof of failure of the test 
than before. 



2 The perturbative framework 

All the methods of the Painleve test [[J, whether differential (Painleve) or al- 
gebraic (Kowalevskaya, Bureau, Conte et al), are based on the following lemma of 
Painleve [§]. 

Lemma. Consider an algebraic differential equation of order N 

E(x,u,e)=0, (1) 

(the boldface characters mean multicomponent) depending analytically on a small 
complex parameter e, defined in the canonical form of Cauchy 

^ = Kb, u, e], x G C, u G C N , eeC, (2) 
ax 

with K assumed holomorphic at e = 0. If the general solution of (|2|) has no movable 
critical points for every e except maybe for e = 0, then 

• the point e = is no exception, i.e. the general solution has also no movable 
critical points there, 
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every coefficient of the Taylor series of u 



+00 

u = V eV") (3) 



n=0 



has no movable critical points. 

A detailed proof of this lemma has been established by Bureau [0, ||, in place of 
the too short proof by Painleve, who merely considered it as an immediate conse- 
quence of the classical theorem of perturbations of Poincare. 

Let us introduce the notation 

+00 

E(x,u,e) = ^£ n EW(i,u(°)...,uW) =0 (4) 

n=0 

in which the equation E^°\x, u^) = is nonlinear and every equation 
E^ n '(x, u.(°\ . . . , u( n )) = 0, n > 1 is linear inhomogeneous in u^. 

Consider an equation E(x, u) = independent of e, a case in which the lemma 
still applies (ref. [|1(| vol. III). The equations E(x, u) = and E^(x, u^ )) = are 



then the same, the equation E^ = is the linearized of E at and equations 
E^ = 0,n > 2, only differ from E^ = in the r.h.s., 

n = 0: E<°> = E(x,u(°>) = 0, (5) 
n = l: E (1) ee E'(x,u (0) )u (1) = 0, (6) 
Vn>2: E (n) ee E'{x, u (0) )u (n) + R (n) (^, u(0) > • • • > u(n ~ 1} ) = °- ( 7 ) 

Suppose one knows a solution u^°) which is global and without movable critical 
points, but which depends on an insufficient number M < N of arbitrary parameters. 
We require, first, that its perturbation (^) in the neighborhood of e = represents 
the general solution, secondly, that each u (n) , in particular in this Letter, be 
without movable critical points. Painleve showed (ref. p. 209 note 1) that M 
of the N solutions of the linearized equation (^) are already known, these are the 
derivatives of u^ ^ with respect to its M parameters, evidently without movable 
critical points. 

Satisfying the first point is not so evident. Indeed, a peculiarity of nonlinear 
differential equations is that they sometimes possess, in addition to a general so- 
lution and its particular solutions, other solutions called singular solutions. Those 
are impossible to obtain from the general solution by giving special values to the 



arbitrary integration constants. For instance, the equation ] 1 



u'" = 2u'u" + 2iu'Wu" -u' 2 -l (8) 
has for general solution: u = e 2ciX+C2 + (cf — l)x/(4ci) +C3 and for singular solution: 



u = — Logcos(a; + C\) + C2, showing the absence of correlation JTTJ] between the 
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structure of singularities of the general solution and that of the possible singular 
solutions. 

Singular solutions are for sure excluded, as requested by the definition of the PP 
and the application of the above lemma, if the equation = is exactly of order 
N. Indeed, the singular solutions, which generalize the notion of envelope, satisfy a 
differential equation with an order smaller than N. 

But the mean used by the available algebraic methods to satisfy this requirement 
is to ask that the equation (D) possesses exactly N Fuchs indices [ for the basic 
vocabulary (singular regular, singular irregular, Fuchsian, non-Fuchsian, essential 
singularity, Fuchs indices) and methods concerning linear differential equations in 
the complex plane, the reader can refer to the classical book of Ince |12|1, chapters 



XV to XVII ] at the movable poles Xq of u^ ) located at a finite distance, which 
implicitly discards possible non-Fuchsian solutions, even if they have no movable 
critical points. 

Let us briefly recall the differences between Fuchsian and non-Fuchsian for our 
purpose which is the question of local singlevaluedness. 

Near a regular singular point x = 0, there exist N linearly independent solutions 

mi 

x A ^^(x)(Logxy, i = l,N (9) 

3=0 

in which the Aj's are complex numbers (the Fuchs indices), rrii positive integers (their 
multiplicity), (pij converging Laurent series of x with finite principal parts. The 
necessary and sufficient condition of local singlevaluedness of the general solution of 
the linear equation is: Aj all integer, no Log terms. 

Near an irregular singular point x = 0, there exist N linearly independent solu- 
tions 

mi 

e Qi ^x Si J2^ij(zi)(^ogx) j , Zi = x 1/q \ i = l,N (10) 

3=0 

in which is a positive integer, Qi a polynomial, a complex number, ip^ a formal 
Laurent series with finite principal part. The question of local singlevaluedness of 
the general solution cannot be settled so easily, because formal series are generically 
divergent. 

In this Letter, we go back to the lemma of Painleve and we implement two 
features. 

The first feature is to also accept that Xq be non-Fuchsian for the linearized 
equation (|J): we require the existence of a fundamental set of solutions of 
EW = which are locally single valued near \ — x ~ x a — 0, and the same property 
for a particular solution of each E^ n ^ = 0, n > 2. 

The second feature needs an additional assumption, namely that the given so- 
lution u(°) be known globally (in closed form). The singular points of equation (|6]) 
can be classified into three types: 

1. the fixed singularities of equation fl5|) located at a finite distance, 
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2. the movable poles xq,xi, ... of u^ ) located at a finite distance, 

3. the point at infinity, which is fixed. 

All three types can be studied because u^ ) is known globally. The first type 
of singularities must not be studied because the PP allows any behaviour around 
them. Each point x ,xi,... of the second type must be studied, with the same 
requirements than in the first feature. As to the third type, namely the point at 
infinity, it should generally not be tested for it is fixed. However, in case there are no 
critical singularities of the first type (fixed critical singularities at a finite distance), 
one can also require the existence of a fundamental set of solutions locally single 
valued near oo; in particular, it may happen the favourable event that, while the 
point xq is non-Fuchsian, the point oo is Fuchsian and allows to conclude about 
global singlevaluedness. 

Our addition to the Painleve test can be detailed as follows. 

1. Assume a given solution u^ ) which is global and has at least one movable pole 
at a finite distance denoted Xq. 

2. Require the linearized equation at u^ ) to have exactly order N, so as to discard 
singular solutions. 

3. Near every movable pole xq of u^ ^ located at a finite distance, build a funda- 
mental set of solutions and require it to be locally single valued. 

4. In case of at least one non-Fuchsian point among xq,Xi, ... at the preceding 
step, and if there is no fixed critical singularity at a finite distance, build a 
fundamental set of solutions near \ = 00 an d require it to be locally single 
valued. 

5. At each higher perturbation order n > 2, build similarly particular solutions 
ir n ) and require the same properties. 

We will call a family Fuchsian or non-Fuchsian according to the nature (Fuchsian 
or non-Fuchsian) of the singular point \ = of the linearized equation. 

The most likely event to occur in our extension of the test, leading to its failure, 
is the detection of a movable logarithm coming from two integer Fuchs indices at 
the point x = oo, of course under the condition of absence of any fixed critical 
singularity at a finite distance. 

The information (pass or fail) which we propose to extract is 

1. not accessible to the method of Kowalevskaya, which accepts the Fuchsian fam- 
ilies as well as the non-Fuchsian ones but does not consider the perturbation 
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2. not accessible to the method of Bureau, which rejects the non-Fuchsian families 
(because the representation by divergent series forbids to conclude) and does 
not consider the perturbation (|3|), 

3. accessible to the method of Conte, Fordy and Pickering, which rejects the 
non-Fuchsian families for the same reason than Bureau and extracts the infor- 
mation from Fuchsian families only at the expense of a Fuchsian perturbative 
computation which may sometimes be long. 

3 An explanatory example: Chazy's class III (N = 

3,M = 2) 

The simplified (i.e. scaled) equation of class III of Chazy [|ll[ 



E{x, u) = u" - 2uu" + 3m' 2 = 0. (11) 

admits the global two-parameter solution 

u — C X~ 2 ~ 6x~\ X — x ~ x 0i (^0) c ) arbitrary. (12) 

For that equation, this solution arises from the local search for all the families of 
movable singularities u ~ Uq\ p \ E ~ E x 9 , X — x ~ x o "~ > represented by Laurent 
series with a finite principal part 

p = -l, g=-4, = -Q X '\ (13) 
p = —2, q = —6, V* = ex 2 — 6% 1 , c arbitrary. (14) 



The linearized equation at (jT3) 

^(ar, u (0) )u (1) = [dl - 2m (0) <9 2 + &ug>d x - 2u^]u m = 0, (15) 



has effectively order N (which means that solution flT2] ) is a particular one, not a 
singular one), it possesses the two single valued global solutions d Xo u^°> and d c u^°\ 
i.e. = x~ 3 , X -2 ) an d ^ has only two singular points x — 0, oo. The point x = is 
irregular singular of rank two for the non-Fuchsian family (|14]) (and regular singular 
for the Fuchsian family (0), with Fuchs indices —4, —3, —2), and it defines a third 
formal solution admitting an essential singularity at x = (ref. |l2j chap. XVII) 



X ^0, c^0: u« = e a /*xMx), (16) 

in which (a, s) denote constants and w(x) a formal Taylor series generically asymp- 
totic, i.e. with a null radius of convergence. This is not the case here, where com- 
putation yields a = —2c, s = —2,w = 1, thus defining the fundamental set of global 
solutions 

V X Vc: u« = X~\ X~\ (e- 2c ^-l + 2c X - 1 )x- 2 /(2c 2 ) (17) 
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and solving the question for the perturbation order n = 1. 

The point x — 00 is therefore, in this too simple an example, useless to study. 
This point is regular singular with Fuchs indices 2, 3, 4, and without any more com- 
putation one concludes to the local singlevaluedness since index 4 cannot generate a 
logarithm. 

Remarks. 

• Going on with the formalism of Painleve's lemma at higher orders constitutes 
the rigorous mathematical framework of the local representation of the general 



solution obtained by Joshi and Kruskal jjTJ 



u = -6 x ~ 1 + cx' 2 (l + z-z 2 /8 + z 3 /lU-7z 4 / 13824 + 0(e 5 )), z = (e/c)e- 2c/x . 

(18) 

This representation reduces to the one given by Chazy (Taylor series in 1/x) 
if one starts from the Fuchsian family 



• The lowering by M = 2 units of the order of the linearized equation is obtained 
by the change of function 

u (i) = x -z v . (q x + 3x -i _ 2cx~ 2 )v" = 0, (19) 

which yields the third global solution by three quadratures. 

The simplified equation ( p]) , which possesses the PP |lT] and therefore for which 
no u^ 1 ' is multivalued, only shows the method. 

We now illustrate on other examples the interest of non-Fuchsian families to 
detect the presence of a movable critical singularity, very often as soon as the first 
perturbation order. 

4 An equation of Bureau (N = 4, M = 2) 

The equation (ref. [Tj| p. 79) 

E(x, u) = u"" + 3uu" - Au' 2 = (20) 
possesses the global two-parameter solution [|J 

u (0) = cx -3 _ 6 Q x -2 ; x = x _ (21) 

for which the linearized equation 

= E'(x, = [dt + 3u^d 2 x - 8u^d x + 3uW]uW = 0, (22) 

whose only two singular points are x = an d X = 00 ? admits the two global 
single valued solutions d Xo u^ and d c u^\ i.e. = x~ 4 , X~ 3 ; leaving only two other 
solutions to examine. 
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For c = it has four global single valued solutions = x~ 5 , X -4 ) X -3 ? X 18 > an d 
from now on we assume c 7^ 0. 

The point x = is singular irregular with rank one, and the two other solutions 
are non-Fuchsian and formally given as 



X 



0, c + : M« = ^V^T^il + O(VX)), (23) 



detecting the presence in (^) of an essential singularity at x — 0, but the generically 
null radius of convergence of the formal series forbids to immediately conclude to 
the multivaluedness of 

The point x — 00 is singular regular, with Fuchs indices —18, 3, 4, 5. The highest 
index, 5, yields a converging series = (l/x) 5 u) 5 (l/x), locally single valued. As 
for the test for the existence of a series = (1/x)~ 18 m>_i 8 (1/x), it detects a 
logarithm arising from the index 3 when one tries to solve the recursion relation for 
the twenty-second coefficient of the series m>_i 8 . 

One thus concludes here as soon as order one, to be compared with order seven 
necessary to ref. ||, after a computation practically untractable without a computer. 

Remark. The lowering by M = 2 units of the order of the linearized equation 
22[) is obtained with 



u« = X - 4 V : (dl - IQx^dx + 3cx~ 3 - 60x~ V = 0, (24) 
and it yields the two other solutions in global form 

c^0:v'! = xAi ? i(24;-3c/x)=X 17/ V 23 (V / 12^), (25) 



v" 



X 17/2 N 23 (Jl2^c), (26) 



where the hypergeometric fonction Fx(24; — 3c/ x) is single valued and possesses 
an isolated essential singularity at x = 0, while the fonction A^ 23 of Neumann is 
multivalued because of a Log x term. 

Remark. The value n = 7 of the perturbation order in e needed by ref. |§ is 
the root of the linear equation n{i m \ n — p) + (i max — p) — —1, with p = — 2,i min = 
— 5,imax = 18, linking the pole order p in the Fuchsian case c = 0, the smallest 
and the greatest Fuchs indices. It expresses the condition for the first occurence of 
a power leading by integration to a logarithm, in the r.h.s. R^ n ' of the linear 
inhomogeneous equation (|7|), r.h.s. created by the nonlinear terms 3mm" — 4m' 2 . 

5 An example in cosmology: Bianchi IX (N = 

6, M = 4) 



The Bianchi IX cosmological model in vacuum [14| is ruled by the sixth order 
system 

a 2 {LogA)" = A 2 -{B- C) 2 and cyclically, a 2 = ±1, (27) 



S 



and it does not possess the PP [T5|, [16| . Let us prove it with our method. 

In the neighborhood of the global solution depending on the four arbitrary pa- 
rameters (ki, k 2 , Ti, r 2 ) [|T7| (x is here denoted r) 

A (o) = j- fc i B (o) = C {o) = j klsmhhjT-n) ^ 
sinh fci(r — Ti) ' sinh 2 A; 2 (r — r 2 ) ' 

the perturbation 

A = A (0) (1 + + 0{e 2 )) and cyclically (29) 

generates a linearized system whose order is indeed equal to iV = 6 (which proves 
that solution ( p8|) is a particular solution, not a singular solution), and the lowering 
by M = 4 units of its order is obtained by the change of function dictated by the 
symmetry of the system: P« = + M« = B<$ - C (1) M 



a 2 A (iy> _ 2A (o) 2 A (i) = 0) (30) 
a 2p(i)'' _ 2 A®B®pW = A(A^B^ - AW 2 )AW, (31) 
a 2 M (D" + 2 (i(o) B (o) _ 2j B(°) 2 )M« = 0. (32) 

Indeed, the four single valued global solutions 

(A®, P (1) ) = d c (Log A<°\ Log(B^ + C^)), c = h, k 2 , n, r 2 , (33) 

are those of the equations (|3~0l)-(j31~D, 

((r-n)cothh(r-ri)- l/h,0), 
(A«,pW + 2^)) = { (0,(r- r 2 )cothfc 2 (r-r 2 ) - 

1 (COth Ki^T — Ti), 0), 

(0, coth k 2 (r - t 2 )), 



and there only remains to study the equation (|32|) , the singular points of which 
(modulo the period of sinh) are r — r 2 = and r = oo. 

At r — r 2 = 0, the equation (|32|) genericaly possesses irregular singular points of 
rank two since the coefficient B^ has there quadruple poles. Its two non-Fuchsian 
solutions are formally (ref. |12| chap. XVII) 



-oo 



T - T2 _> : M« = e Q /^ £ A fc (r - r 2 ) fe+s , A ^ 0, (35) 



fc=0 



with 



a: 



±.2k 1 1 sinh A;i(r 2 — Ti), s = 1 =F 2 cosh A;i(r 2 — Ti). (36) 



The two generically irrational values for the exponents s allow to conclude only if the 
Taylor series Afe(r — r 2 ) k can be summed, which is the case at least for k\ = k 2 = 
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where the two solutions are globally known 



h = k 2 



: 




M (1) = 0, t 



z 4 e 4z dz, 



r-r 2 
T1-T2 



(37) 



The second solution implies the presence of a logarithmic branch point at t = 0, or 
at t = 00 as well. This singularity persists for (ki,k 2 ) 7^ (0,0) and this proves the 
absence of the Painleve property for the Bianchi IX model in vacuum. 

6 Conclusion 

The main application of this extension to the Painleve test is the case of a 
non-Fuchsian family: if the unperturbed solution is known in closed form and if 
the linearized equation has no fixed singularity at a finite distance, the study of the 
point at infinity often allows to conclude to non-integrability (in the sense of absence 
of the PP) much more rapidly than the above mentioned algebraic methods. The 
present method constitutes an algorithmic extension to the Painleve test. 

For the two examples given where the Painleve test fails, this failure was already 
known, and we do not know yet if there exist equations where a failure would only 
be detectable by the present extension and, of course, by the a— method. 

Nevertheless, the algorithmic cost of the present method is much smaller than the 
cost of earlier methods, and this is particularly sensitive on the example of Bureau. 

Moreover, in our two examples, the perturbation order n — 1 is sufficient to 
dectect a failure, which the perturbative Fuchsian test detects at the respective 
orders n = 7 (Bureau) and n = 2 (Bianchi IX). Despite the examination of several 
other equations without the PP and admitting a non-Fuchsian family, we could not 
find an example requiring an order n greater than one to exhibit a failure of the 
test. Such a feature may be generic and certainly deserves future investigations. 

Acknowledgements. We thank D. Bessis, A. Latifi, A. Mezincescu and A. Pick- 
ering for fruitful discussions. 
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